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The critical dynamics of superconductors in the charged regime is reconsidered within field-theory. For the 
dynamics the Ginzburg-Landau model with complex order parameter coupled to the gauge field suggested 
earlier [Lannert et al. Phys. Rev. Lett. 92, 097004 (2004)] is used. Assuming relaxational dynamics for both 
quantities the renormalization group functions within one loop approximation are recalculated for different 
choices of the gauge. A gauge independent result for the divergence of the measurable electric conductivity 
is obtained only at the weak scaling fixed point unstable in one loop order where the time scales of the order 
parameter and the gauge field are different. 
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1. Introduction 

The nature of the static phase transition in superconductors was an open question for decades, 
since due to the large correlation length of existing superconducting materials the effect of critical 
fluctuations was hard to observe in the vicinity of the critical temperature T c . Moreover it was 
unclear if the phase transition is of first or second order and, if it is of second order, to which 
universality class the transition belongs. The appearance of high-T c superconducting compounds 
with short correlation lengths [ [1] made the critical region of superconductors experimentally 
accessible. In turn this leads to comparison of the experiment with analytical results in order to 
establish the critical properties of superconductors and their universality class. 

From the theoretical point of view the static critical properties are now well understood. The 
theoretical model for the description of the static critical properties was formulated in [ [2] and 
contains, besides the two component (n = 2) order parameter (OP) and its fourth order interaction 
term, a minimal coupling to a gauge field (GF) due to the charged character of the OP (macroscopic 
wave function of the Cooper pair). The coupling of the OP to the GF introduces an essential 
difference to the superfluid phase transition, where the OP is uncharged. 

While for type-I superconductors fluctuation effects are weak and a mean-field analysis can be 
applied, the situation with type-II superconductors is more complicated, since here fluctuations of 
the OP cannot be neglected. The first renormalization group (RG) analysis in one-loop approxi- 
mation [[2] lead to the conclusion that a stable fixed point (FP) and therefore a second order phase 
transition only exists for OP dimensions n larger than 365.9. For OP dimension n = 2 no stable 
FP was found and the runaway solution was interpreted as a weak first order phase transition. A 
two loop order calculation of the renormalization group equations within field theory indicated the 
possibility of a continuous phase transition for n = 2 [ [3] if certain resummations are performed. 
This has been investigated further and led to the following picture [S]: There are four FPs, two 
uncharged ones known from the standard Ginzburg-Landau- Wilson (GLW) model and two charged 
ones. Depending on the initial (background) conditions for the flow of the fourth order coupling 
and the coupling to the charge (defining the value of the Ginzburg parameter n) one obtains a 
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runaway flow or the flow reaches the stable charged FP. The other charged FP is reached starting 
on the separatrix (defined by k = k c ) separating the attraction region of the charged FP from the 
runaway region. Physically this FP reached on the separatrix describes tricritical behavior indi- 
cating that a tricritical point separates the first order transition of the superconductors of type I 
from the second order transition for superconductors of type II. Results of duality arguments [ [7] 
and Monte Carlo calculations [ 8] are in agreement with this picture, which is also supported by 
experiment [l9| I10j. 

An important question in calculating critical properties like the critical exponents is their 
dependence on the gauge used for the GF. Physically observable quantities should be independent 
of the specific gauge used. This has been shown in one loop order [[5]. However non observable 
quantities might depend on the gauge. Thus the critical exponents v or a of the penetration length 
A (proportional to the correlation length £ of the OP correlation function at the charged FP) and 
the specific heat respectively should be independent of gauge, whereas the critical exponent of the 
static OP correlation function turns out to be gauge dependent [[6]. 

The issue of dynamical critical properties was less studied. Experimental investigations give no 
consistent picture of the dynamical critical exponent z (for references see Ref. [E]) obtained from 
measurements of the electrical conductivity. The values of the dynamical critical exponents found 
vary between 1.5 and 2.3. 

Theoretical predictions of z are mainly based on using the results of known universality classes 
(model A, C, E or model F, for a review see [Q2])- Model A is the simplest model assuming a 
relaxation equation for the OP without coupling to other conserved densities, whereas model C 
contains a coupling to a diffusion equation for a conserved density (only relevant if the specific heat 
is diverging) . Models E and F define the universality class of the critical dynamics of the superfluid 
transition in 4 He. Dynamic equations from the vortex-loop model [ [131 [14] indicate a relation to 
that universality class. So far no systematic derivation of dynamical equations with mode coupling 
terms derived from Poisson bracket relations using for instance methods described in [ 15J has been 
performed. 

Monte Carlo simulations [ [13] in the limiting cases of large and small values of the Ginzburg 
parameter found values of the dynamical exponent of 2.7 and 1.5 respectively. Using the vortex 
model [[TT] for the superconducting transition the dynamic critical exponent has been calculated 
analytically in the two limits mentioned above. The resulting exponents were 5/2 and 3/2 respec- 
tively and they were related to the dynamic exponents obtained in simulations. A discussion [ [16] 
arised on the reason of the deviation of the Monte Carlo results from the expected value of model 
A, which would lead to a dynamical exponent of roughly 2. Recently in Ref. [[14] for extreme type 
II superconductors, where the uncharged FP describes the static behavior, it was argued using 
scaling and duality arguments that z = 3/2 exactly. 

A concrete dynamical model defined by equations of motion for the OP and the GF has been 
presented in Ref. [ [T7] . Both equations are of relaxational character and the essential parameter is 
the ratio w — /T a of the kinetic coefficients of the OP and Ta of the GF. A dynamic FP, w* 
finite, implying strong scaling with a common dynamical critical exponent z for the OP and the 
GF in the Feynman gauge (one adds a quadratic term in the divergence of the GF to the static 
functional) has been found in one loop order. However the question of the gauge dependence of 
the dynamical exponents was not addressed in Ref. [117], 

Another important issue not considered so far in the discussion of the dynamical critical behav- 
ior concerns the time scales entering the problem. One has to discriminate between the time scale 
for the OP and its dynamic critical exponent z^ and the time scale of the GF and its dynamic 
critical exponent za- They are defined by the characteristic frequencies of the dynamic correlation 
functions for the OP and the GF respectively 

city ~ k^ gi ,{k£) uja ~ k ZA g A (kO . (1) 

and <ju are scaling functions which are finite and nonzero at the superfluid transition T c . 
The dynamic critical exponents are different [ [18] if weak scaling holds otherwise one says 
the dynamic correlation functions obey strong scaling. So far in the discussion the last case was 
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assumed. Even if for the uncharged FP the dynamics would be in the universality class of model 
E, it is the weak scaling FP which is stable [I19j. 

Only one of these dynamical exponents is observable namely za entering the frequency depen- 
dent diverging electric conductivity 

<r(£,fc = 0,w) ~ ( 2 ) 

r/A = 4 — d is the anomalous dimension of the GF independent of the gauge [ 20j . This implies that 
the dynamical critical exponent Za has to be gauge independent. No such condition holds for z^. 
If therefore it turns out that z^ is gauge dependent strong scaling cannot hold. 

In this paper we apply the field-theoretical RG approach in the minimal subtraction scheme and 
dimensional regularization to the dynamical model of Ref. [ [17] in order to investigate the gauge 
dependence of the critical dynamics. It turns out that the dynamical exponent z$ of the OP is gauge 
dependent and a gauge independent value for za is only obtained in the weak scaling FP w* — > oo, 
which however has been found unstable. Therefore one has to question the physical relevance of 
the one loop result. This goes together with the problems in statics of the Ginzburg-Landau model 
at the one loop order level. 



2. Model 

Static critical properties of charged superconductors in d dimensional space are described by 
the Abelian Higgs model with the static functional [G]: 

r i i n/2 ° 1 

2—1 

with the complex n-component OP ipo (generalized to n/2-components, the superconductor being 
the case n = 2) and the rf-dimensional GF A . The bare parameter r is proportional to the distance 
from T c and the bare coupling e is the effective charge, u is the usual fourth order OP coupling. 
In all subsequent calculations we add to the static functional the contribution 

H = ~(V-A ), (4) 

which allows to choose the gauge. In the limit <f — * the Landau gauge is recovered, while for 
system with ? — > 1 the Feynman gauge is achieved. 

The OP and the GF are nonconserved quantities thus their dynamical behavior may be de- 
scribed by relaxation equations. Such a set of equations of motion has been suggested [[17]: 

fofoM _ 9 f SH . a . 
dA . a . 5H , „ 

-1 A" 



dt SA 0tOl 

The OP relaxes with the kinetic coefficient while the GF relaxes with the kinetic coefficient 
The stochastic forces in © are related to the kinetic coefficients and satisfy the relations: 

<0 i (x,t)0+(x / ,t / ) > = 4f v ,S(x-x')r5(*-i')<% ! 

<6i(x,t)> = 0, (6) 

<6 a (x,t)0p(x!,lf)> = 2fA5(x- X ')5(t-t')5 a0 , 

<e a (x,t)> = o. 
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The indices i,j adopt the values 1, . . . , n/2, while the indices a, (3 adopt the values l,...,d. 

As it was pointed out in Ref. [[17] the equation for the GF in J5]) can be derived from Maxwell's 
equations in their low-frequency form. In this case the inverse transverse coefficient T^ 1 for the 
GF can be identified with the bare normal electrical conductivity. 

We study this critical dynamics by applying the Bausch-Janssen- Wagner approach [ [21] of 
dynamical field-theoretical RG. In this approach, the critical behavior is considered on the basis 
of long-distance and long-time properties of the Lagrangian incorporating the features defined by 
the dynamical equations of the model. In order to keep the powers in the interaction terms low 

an auxiliary OP density V> and an auxiliary GF A have been introduced. The dynamic equations 
defined by ©-([6]) are then described by the unrenormalized Lagrangian £ = Cq + L\ separated 
into a Gaussian part 



£ ■ = / d d xdt [ - 4f ^ #m$<M ~ 2 ^a K ( — + tj, (r-A)j Vo,i + 

"'—i i— 1 



n/2 

E^o,i ( ^+f>(f-A)j 00+ + Ac ■ + f A (V(1 - J)(V • Ao) - AA ) ) , (7) 



i=i 



and an interaction part 

n/2 n/2 



z— 1 i— 1 

n/2 ^ n/2 



i=l i=l 

n/2 

f ^e 2 A 2 J^o^ + + f Ae 2 A • A |</f ] . (8) 

z=l 

3. Perturbative expansion and renormalization 
3.1. Vertex functions 

In order to proceed the vertex functions for the OP and the GF are calculated in one loop 
expansion. We keep the general structure of the two point vertex functions [[32] and separate the 
static contributions. The calculation leads to the following general form of the OP vertex function 

= -twft^ + +2TVf& + . (9) 
In one loop order the dynamical OP function reads 

fU + = 1 • 1< J I , /' — : (lr {kk ' ? 



W+ (r+(k+k') 2 )k' 2 (-iuj+f^(r+(k+k') 2 )+f A k' 2 ) V fc ' 2 

( -. . |, (1 Q) 

whereas the static OP vertex function is given by 



2 \ T ~ ^ 6 J (f + k»y y v 7 /t' 2 

.ay 1 A 2 (^') 2 1 ,((2fc + fc')fc') 2 



46 ' (r + (jfc + k') 2 )k' 2 \ k k' 2 + ^ k' 2 ' (1L) 
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The same general structure holds for the GF vertex function 



ff A = -i^f A + r A t^, (12) 

with the one loop expression 

fr* = r> + t^[ (k + 2kT(k + 2kr 

2 J (f + k' 2 ){f + (fc+fc') 2 )h^ + r V; (r + k' 2 ) +T 4 ,{r+ (fc+fc') 2 )] 
for the dynamic GF function and 

k a kP 1 fc a £^ 



= k 2 [S a(3 - 



k 2 <; k 2 

09 f 1 /c 1 (fc + 2fc') Q (fc + 2fc') /3 \ , , 

for the GF static vertex function. 

3.2. Renormalization and field theoretic functions 

In order to get finite results we perform the renormalization within the minimal subtraction 
scheme introducing renormalization factors leading to the renormalized parameters. The renormal- 
ization factors of the GLW part of the static functional are introduced as usual 

ipo,i = Z^/ 2 ipi, r - r c = Z^ 1 Z r r , u = k e A^ 1 Z u Z^, 2 u , (15) 

where the shift of the phase transition temperature f c has been taken into account, n represents 
the free wave vector scale and e = 4 — d. We also introduced the usual geometric factor 

A d = r(i-i)r(i ■ e ^ nd 



2/ (2ir) d 



2/ V 2/ (2tt) 

where fid is the surface of the ri-dimensional unit sphere. The additional renormalization factors 
due to the presence of the GF and its coupling to the OP are introduced as 

A , a = Z A /2 A a , e 2 = n £ A d : 1 Z e 2Z- 1 Z A 1 e 2 , $ = Z- l Z A <;. (16) 

From Ward identities one derives the relations [\E\ 

Z e2 =Z^, Z, = l. (17) 

which show that only one additional (for the GF) renormalization constant with respect to the 
GLW model appears. 

In dynamics two additional renormalization factors are needed which are 

^ = Zl'%, A Q . a = Z x [ 2 A a . (18) 
For the kinetic coefficients no new factors are necessary 

r</> = = z^ 2 z~ 12 Ttp, (19) 

t A = Z Fa T a = Z\ /2 Z A 1/2 T A , (20) 

where the second relations are due to the structure of the vertex functions (see ([9]) and lfl2|) ). 
The renormalization factors calculated in one loop order for the statics are 

2 2 

^ = l + (3-?) — . Z A = \-n e —, (21) 

£ 0£ 



5 



M. Dudka R. Folk G. Moser 



„ n + 2u e 2 .. .. . . 

Z r = l + — c-, Z u = l + — + 2?-, 22 

bee 6 e u e e 

as presented in [ [3] . The dynamic renormalization factors are 

Z i = WS-O-^^^-, (23) 



e e 1 + w 
2e V3 to, 



^ = 1+nj; r - , (24) 



where we have introduced the time scale ratio 

(25) 

whose FP value determines if strong (w* ^ or oo) or weaA; scaling (w* = or oo) holds. 

From the Z-factors one obtains the C-functions leading to the /3-functions which determine the 
FPs. The critical exponents describing the critical properties are then expressed by the values of 
the C-functions at the stable FPs. We use the unified definition 

, r > d In Z„ , 

«< a >''> = -lE7- (26) 

for all ^-functions, where a denotes either any model parameter {a} = {u, e 2 , T^, Ta, w} or any 
density tjj, ip, A, A. The static C-functions following from lj2Tj) and lf22]) read 

= (3-0e 2 , U = ~n^, (27) 
n + 2 

Cr = -g- « - ^ 2 , ^ - (3 - ?)e 2 , (28) 

C = ^u-2,e 2 + l8-. (29) 
6 it 

The C-function of the gauge parameter follows from (fl6| and (fl7| as C? = 0. The dynamic £- 
functions follow from (|23|) . They are 

and the C-functions for the kinetic coefficients are obtained from the above relations lf!9|) and (|20|) 

Cr«, = ^(Cv _ C^)> Cr^ = - Ci)- (31) 



Inserting the one loop results lj30j) they read 



According to Eqs. (|25f and l[3"Tj) . the C-function of the time scale ratio w is then found as 

? n /l 1 
1 + tu + 2 V 3 2^ 



^=Cr,-Cr.=e 2 (3--^- + ;fi- 7 l)) . (33) 



4. Fixed points and exponents 

The behavior of the model parameters under renormalization is described by the flow equations 

^=/WW,0, (34) 
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where the right hand sides of the equations are determined by appropriate /3-functions (the index 
i runs over the set of parameters) . The /3-functions are generally defined as 

({"0,0 = on [-a -piCj, - qi(A + CqJ , (35) 

The coefficients Cj , and & follow from the general renormalization of a, as used above 

&i = K^A^^Z-^Z^Z^ai . (36) 

The FPs {a*} are defined by the zeros of the right hand sides of Eqs. J34|). A FP is stable if all 
stability exponents u>i are positive. The stability exponents oj, are defined by the eigenvalues of the 
matrix d/3 ai /daj\ { a y={ a *}. 

4.1. Statics 

The one-loop static /3-functions appear to be independent from the gauge determined by the 
parameter c; and are the same as in Ref. [[2]. 

P u = f-eu - 6e 2 u + ^-^u 2 + 18e 4 ^ , (3 e 2 = e 2 (-£ + ^e 2 ) . (37) 

These /3-functions have four FPs: (i) the Gaussian FP with u* = e* 2 = 0, (ii) the FP of the 
uncharged XY- model (it* = 6e/(n + 8),e* = 0), (iii) the tricritical FP with u* = 0, e* 2 — 6e/n 
and (iv) the charged FP with both couplings u* ^ and e* 2 = 6e/n nonzero. In one loop order 
the charged FP exists only for n > 365.9 [|2|. However in higher loop order and using summation 
procedures a charged FP is also found for n = 2 [E2]- In one loop order for n < 365.9 starting 
with e 2 ^ and u ^ the flow escapes to e 2 — * e* 2 and u — > oo. 

Critical exponents are calculated at the stable accessible FP. Expressions for the correlation 
length critical exponent v , the specific heat exponent a and the pair correlation function critical 



exponent rj are as follows: 

n + 2 

y~ = 2-Cr(u*,e*,0 + Ctf(«*,eV) = 2 —u* + 3e* 2 ; (38) 

a = 2~dv; (39) 
7] = C^K,e*,0-(3-0e* 2 - (40) 



All other static critical exponents may be found from scaling relations. Since the FP values u* 
and e* are gauge independent this is also valid for v and a although the gauge ? appears in the 
(■-functions Eq. ((3Sj) explicitly. In contrast this dependence on the gauge remains in rj (see Eq. 
(|40j) ). All the other exponents become gauge dependent since they are in any case related to r] 
by the scaling laws. However note that only at the charged FP (e* 2 ^ 0) a gauge dependence 
enters the (-functions. This holds in all orders of the loop expansion due to the structure of the 
interaction part of the Lagrangian ([8]) . 

From the renormalization of <r in Equ. lf!6|) the flow of the gauge parameter reads 

(«) 

reaching zero at the charged FP. This means only the transverse gauge is invariant under renor- 
malization. 

4.2. Dynamics 

The FP values for the time scale ratio can be found from the zeros of the corresponding (3- 
function 

A. = «(fr, - Cr J = e 2 W ( 3 - , + ^ + \ (| - 2.)) . (42) 
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These zeros depend in one loop order on the value of the minimal coupling (charge) e at the stable 
static FP. Only for non-zero charge a specific value for the possible FP is found in this order (w 
has to be positive) 



t n - 12(3 - ?) + J(n - 12(3 - <r)) 2 + 24n(18 + n) 

w = — — - (43) 



4(18 + n) 

The corresponding stability exponent oj w 

9 P „*2„..* 



" dw ~ e * w \{\ + w*y + 4w* 2 ) (44) 

turns out to be positive and thus the finite FP is stable. Another FP is the infinite FP. In order 
to study its stability we introduce 

'-IT= < 45 » 

mapping the range of w into the interval [0, 1]. The corresponding /3-function reads 
0, = P(l - P)e 2 (3 - ? + CP + \ (~ - ^) ) = (1 - P)e 2 ((3 - ,)P + + - J) ■ (46) 

The FP p* = 1 corresponding to the infinite FP w* -1 = however is unstable with the stability 
exponent uj p = —(1 + 18/n)e. 

The dynamical critical exponents for the OP and the GF are calculated from 

Z0 = 2 + Cr*(u*,e* J u;*,O, z A = 2 + Cr A (u*,e*,w*) . (47) 

Whereas an explicit gauge dependence is found in (r^ this is not the case for Cr A at least in one 
loop order. It is expected that this property holds in all orders. At the charged FP (e* 2 = 6e/n) 
in general the exponents take the values 

„ 18 6 e 
Z * = 2+ - £ - ? nTT^ ; (48) 

ma = (49) 

Since w* depends on the gauge for the finite FP, the dynamical exponent of the GF, Za, depends 
on the gauge in contradiction to the physical requirement that observable quantities should be gauge 
independent. Only in the case of the infinite FP a gauge independent value is possible. Then the 
OP exponent would be finite 

18 , s 

= 2 + — e, (50) 

and also independent of the gauge but different from z A - However the infinite FP is not stable in 
one loop order. It should be noted that for the gauge ? = 1 the results of [[17] are reproduced. It 
should also be remarked that in one loop order the calculation of the field theoretic functions is 
simpler in the transverse gauge (c = 0) due to the observation that less graphical contributions in 
the loop expansion are nonzero [[5j. The gauge independence of the static FP values for the minimal 
coupling and fourth order coupling has not been proven to our knowledge. The finite FP value of 
the dynamical time scale w is already gauge dependent in one loop order, as has been shown here, 
leading to a gauge dependent critical dynamic exponent of the GF. This in our opinion would 
hold for every order for the finite nonzero FP value w*. At the uncharged FP both dynamical 
exponents would be independent of the gauge. 
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5. Conclusion and Outlook 



We have demonstrated for the dynamical model suggested by Lannert, Vishveshwara and Fisher 
[ [17] that the OP dynamical exponent is gauge dependent. Therefore also the exponent for the 
divergence of the electric conductivity is gauge dependent at the stable strong scaling FP where 
the dynamical critical exponents of the OP, z^, and the GF, za, are the same. A way out in a 
higher loop order calculation is only possible when the stability of the strong scaling FP is changed 
and the infinite weak scaling FP becomes the stable one. As one knows from other examples one 
loop order calculations may not be conclusive and lead to results which have to be taken with 
care. In higher order perturbation expansion the stability of the FPs may be changed. A two 
loop calculation may clarify the situation. Recent progress [[24] in the nonperturbative version of 
dynamic renormalization theory applied to this more complicated model would be also worthwhile. 
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